Search for pulsed 7-ray emission from globular cluster M28 
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ABSTRACT 

Using the data from the Large Area Telescope on board the Fermi Gamma-ray Space Tele- 
scope, we have searched for the 7-ray pulsations from the direction of globular cluster M28 (NGC 
6626). We report the discovery of a signal with the frequency consistent with that of the ener- 
getic millisecond pulsar (MSP) PSR B1821-24 in M28. A weighted H-test test statisic (TS) of 
28.8 is attained which corresponds to a chance probability of ~ 10 -5 (4.3a detection). With a 
phase-resolved analysis, the pulsed component is found to contribute ~ 25% of the total observed 
7-ray emission from the cluster. On the other hand, the unpulsed level provides a constraint for 
the underlying MSP population and the fundamental plane relations for the scenario of inverse 
Compton scattering. Follow-up timing observations in radio/X-ray are encouraged for further 
investigating this periodic signal candidate. 

Subject headings: gamma rays: stars — pulsars: general — stars: individual (PSR B1821-24, PSR J1824- 
2452A) 



1. INTRODUCTION 

The first millisecond pulsar (MSP), which is 
a rejuvenated old neutron star through accret- 
ing matter from its companion, was discovered 30 
years ago (Backer et al. 1982). It has long been 
suggested that they are the descendants of low- 
mass X-ray binaries (LMXBs; Alpar ct al. 1982). 
In comparison with the Galactic field, the for- 
mation rate per unit mass of LMXBs in globu- 
lar clusters (GCs) is orders of magnitude higher 
because of the frequent stellar encounters (Katz 
1975; Clark 1975; Pooley et al. 2003; Hui, Cheng 
& Taam 2010). Therefore, it is not surprising 
that GCs should host a large population of MSPs. 
Since the first cluster MSP, PSR B1821-24, has 
been discovered in M28 (Lyne et al. 1987), dedi- 
cated radio pulsar surveys towards different clus- 
ters have resulted in the currently known popula- 
tion of 144 MSPs in 28 GCsH 

Since the launch of Fermi Gamma-ray Space 
Telescope, a new population of 7— ray emitting 



5 see |http: //www. naic.edu/ ^pfreirc/GCpsr, html for up- 
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globular clusters (GCs) have been detected (Abdo 
ct al. 2010a; Tarn et al. 2011). As MSPs arc the 
only known steady 7-ray sources in GCs, they are 
suggested to be the contributors for the observed 
emission. The 7-rays from a GC are interpreted 
as the collective contribution from the entire pul- 
sar population resides in it (Abdo et al. 2010a; 
Cheng et al. 2010; Hui ct al. 2011). For most of 
the cluster MSPs, detecting the 7-ray pulsations 
is very challenging as the 7-ray flux of individual 
pulsar is weak. To exacerbate the situation, the 
spatial resolution of Fermi Large Area Telescope 
(LAT) docs not allow individual MSPs in a GC 
to be resolved. This results in a high background 
that make the pulsation search for any individual 
MSP very difficult. However, there are two no- 
table exceptions, namely the GCs NGC 6624 and 
M28. Both of these GCs contain a very energetic 
and young MSP so that they can possibly stand 
out from the background. 

For PSR J1823-3021A in NGC 6624 (P = 
5.44 ms), its large spin-down rate, P = 3.38 x 
10~ 18 s s _1 , implies that it is the youngest MSP 
(r ~ 25 Myrs) ever detected (Freire et al. 2011). 
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Its spin-down luminosity is E = 8.3 x 10 35 erg s 1 , 
which is ~ 1 — 2 orders of magnitude greater than 
the typical MSPs in GCs (Bogdanov et al. 2006). 
Thanks to the accurate timing model provided 
by the dedicated radio observations, its 7— ray 
pulsations has been revealed by Fermi (Freire et 
al. 2011). Through a phase-resolved analysis, it 
has been shown that this single pulsar dominates 
all the observed 7— rays from NGC 6624 (Freire 
et al. 2011). The 7— ray conversion efficiency 
of PSR J1823-3021A, L 7 '/E ~ 0.1, is found to 
be comparable with other 7— ray detected MSPs 
(Abdo et al. 2010b), where L 7 is the 7— ray lumi- 
nosity. 

PSR B1821-24 in M28 (hereafter M28A) is 
very similar to PSR J1823-3021A in many as- 
pects. Its period (P = 3.05 ms) and spin-down 
rate (P = 1.61 x 10~ 18 s s _1 ) imply its age 
and spin-down power to be r ~ 30 Myrs and 
E = 2.2 x 10 36 erg s _1 respectively, which makes 
it as the most energetic MSP has been found so 
far (Bogdanov et al. 2011). Together with its non- 
thermal X-ray spectrum, its sharp and narrow X- 
ray pulse profile strongly indicate that most of the 
observed X-rays from M28A are originated from 
the magnctosphcre. 7— ray emission from M28 has 
been detected by Fermi LAT (Abdo et al. 2010a). 
As its Lj is only a fraction of the spin-down power 
of M28A, it is possible that this pulsar can have 
a significant contribution to the observed 7— rays. 
Together with its relatively short distance, 5.5 kpc 
(cf. Harris 1996; 2010 version), the cluster M28 is 
a promising target for searching 7— ray pulsation. 
In this Letter, we report our recent search for the 
possible pulsation from this GC by using Fermi 
LAT data. 

2. DATA ANALYSIS AND RESULTS 

In this work, we used the Fermi LAT data 
between 2008 August 04 and 2012 January 31. 
For the data analysis, the Fermi Science Tools 
v9r23pl package, available from the Fermi Sci- 
ence Support CentciH, was used. We used Pass 
7 data and selected events in the "Source" class 
(i.e. event class 2) only. In addition, we ex- 
cluded the events with zenith angles larger than 
100° to greatly reduce the contamination by Earth 
albedo gamma-rays. The instrumental response 

http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/software/ 



functions (IRFs) "P7SOURCE.V6" were adopted 
throughout the study. Events were selected within 
a circular region-of-interest (ROI) with a diame- 
ter of 10° centered at the optical center of M28. 
Photon energies are restricted in the range of 
200 MeV-300 GeV. This set of cuts is adopted 
throughout this work. 

To investigate the spectral characteristic of M28 
with the updated IRFs and background model, we 
performed an unburnned likelihood analysis with 
the aid of gtlike by assuming a point source with 
power-law with exponential cutoff (PLE) of the 
form dN/dE oc E~ r exp (— E/E cuto s) at the nom- 
inal position of M28, where T and £^ C utoff are the 
photon index and the cutoff energy respectively. 
For modeling the background, we included the 
Galactic diffuse model (gal_2yearp7v6_v0 . f its), 
the isotropic background (iso_p7v6source.txt), 
as well as all point sources reported in the 2FGL 
catalog within 10° from the center of the ROI. 
All these 2FGL sources were assumed to be point 
sources which have specific spectrum suggested by 
the 2FGL catalog (Nolan et al. 2012). While the 
spectral parameters of the 2FGL sources locate 
within the ROI were set to be free, we kept the pa- 
rameters for those lying outside our adopted ROI 
fixed at the values given in 2FGL (Nolan et al. 
2012). We allowed the normalizations of diffuse 
background components to be free. The best-fit 
PLE model is characterized by T = 0.96 ±0.22 and 
E cutoS = 1.41±0.3 GeV with a test-statistic (TS) 
value of 825 which is highly significant. We have 
tested the robustness of the spectral results by re- 
peating the analysis with different size of ROI. The 
fitted parameters from independent analysis are 
consistent within la uncertainties. In this model 
the photon flux between 200 MeV & 300 GeV 
was found to be (2.39 ± 0.22) x 10" 8 cnT 2 s _1 . 
The corresponding integrated energy flux is / 7 = 
(3.17 ± 0.29) x 10" 11 erg cm" 2 s" 1 . The spectral 
results are consistent with those reported by Abdo 
et al. (2010a) within la uncertainties. 

Assuming M28A is major contributor for the 
7-rays from M28, we search for the possible 7-ray 
pulsation from this GC. We started by adopting 
the timing ephemeris determined for M28A from 
a long term observation with Rossi X-ray Tim- 
ing Explorer (Ray et al. 2008), which are tab- 
ulated in Table 1. For barycentric correction, we 
used the updated planetary ephemeris JPL DE405 
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throughout this analysis. Following the method 
proposed by Kerr (2011), we used the best-fit 
model resulted from the phase-averaged likelihood 
analysis as described above for assigning weight to 
each 7-ray photon by computing the probability 
that it originates from M28. This enables us to 
alleviate the problem of source confusion with a 
more efficient background rejection. We then as- 
sign a pulsar spin phase to every gamma-ray pho- 
tons with energies >0.2 GcV and within 5° from 
M28A's direction (see Tab. 1). A promising signal 
with a weighted H-test TS of 28.8 has been found 
by directly folding up the data with this ephemeris 
(de Jager & Biisching 2010). The folded 7— ray 
pulse profile and the phaseogram (i.e. pulse phase 
as a function of time) with the weighted photons 
are shown in the upper panel and the lower panel 
of Figure Q] respectively. 

According to Figure [TJ it appears to have two 
peaks with one broader than the other. We then 
define the phase-interval for the peak 1 and peak 2 
to be — 0.4 and 0.55 — 0.75 respectively. The rest 
is defined as the off-pulse component. With this 
definition, we show the Fermi LAT count maps of 
the sky region around M28 at different phases in 
Figure [2J During the on-pulsc intervals, a point- 
like 7— ray source can be clearly seen at the pulsar 
position which is illustrated by the yellow cross. 
On the other hand, a faint diffuse excess is found in 
the off-pulse phase (i.e. 0.40< <0.55 and 0.75< 
<f> <1.0). However, the limited photon statistic 
does not allow us to constrain the extent of this 
putative feature. 

This on-and-off nature of the 7— ray emission 
from M28 provides a strong support for the pres- 
ence of periodic signal and leads us to a more de- 
tailed investigation. For investigating the possi- 
ble spectral variations among peak 1, peak 2 and 
the unpulsed component, we performed a phase- 
resolved likelihood analysis. We have fitted their 
spectra with both simple power-law (PL) of the 
form dN/dE oc E~ r and PLE. The results arc 
summarized in Table 2. According to the PL 
fits, there is no obvious change of the spectral 
steepness. We notice that the likelihood analy- 
sis that incorporate the PLE model results in a 
higher TS for all three components. For peak 2, 
we found that the spectral parameters for the PLE 
fit cannot be properly constrained. Therefore, we 
fixed the photon index at the value inferred in the 



phase-averaged analysis (i.e. T = 0.96). Same as 
the cases of the PL fits, within the tolerence of the 
statistical uncertainties, we do not find any con- 
clusive evidence for the spectral variation across 
the phase. Assuming the off-pulse component has 
a constant contribution across the whole phase, 
~ 75% of the total observed flux is originated from 
this component. 

3. SUMMARY & DISCUSSION 

In this Letter, we report our detection of 7— ray 
pulsation from the direction of the GC M28. We 
have found a periodic signal which is presumably 
originated from its energetic MSP M28A. Based 
on our phase-resolved analysis, the pulsed com- 
ponent contributes ~ 25% of the total observed 
7— rays. At a distance of d = 5.5 kpc, this im- 
plies an on-pulse luminosity of L 7 = And 2 /o/ 7 ~ 
3 x 10 34 /si erg s _1 , where fa is the fraction of 
the sky covered by the 7— ray beam. Assuming 
the pulsed emission is originated from M28A, this 
suggests a 7— ray conversion efficiency of L-y/E ~ 
0.01/n. Some of the MSPs in the Galactic field, 
such as PSRs J2124-3358 and J0437-4715, have 
their L 1 /E found at this level (Abdo ct al. 2010b). 
However, this is lower than L 7 /E ~ 0.08 as de- 
rived from the nearby MSPs (Abdo et al. 2009). If 
one adopt this as the intrinsic 7— ray conversion ef- 
ficiency of M28A, this might suggest the observed 
period derivative is largely dominated by the ac- 
celeration of the pulsar along the line-of-sight due 
to the gravitational field of the cluster. 

On the other hand, the off-pulse luminosity is 
found at the level of L 1 ~ 8 x 10 34 erg s" 1 . This 
estimate is useful for constraining the collective 
properties of the rest of the MSP population in 
M28. For explaining the unpulsed 7— ray emis- 
sion from GCs, there are two main theories. One 
stream interprets the 7— ray emission from a GC 
is originated from the collection of the magne- 
tospheric radiation from the entire MSP popula- 
tion resides in it (Abdo et al. 2010a; Venter et 
al. 2008,2009). Assuming an average spin-down 

power of (e^J ~ 2 x 10 34 erg s" 1 and a character- 
istic conversion efficiency of ~ 0.08, the off-pulse 
luminosity enables us to estimate the number of 
the rest MSP population to be Amsp ~ 50. This 
suggests that about one fifth of the underlying 
population has already been uncovered. 
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Apart from the aforementioned standard sce- 
nario, inverse Compton scattering (ICS) between 
the relativistic pulsar wind particles and the ambi- 
ent soft photons has also been proposed as another 
possible explanation for the origin of the 7— ray 
from GCs (Bednarek & Sitarek 2007; Cheng et al. 
2010; Hui et al. 2012). Cheng et al. (2010) found 
that the observed 7-ray spectra of GCs can gen- 
erally be well- modeled by ICS between the e~/e + 
in the pulsar wind of the whole MSP population 
in a GC and the Galactic background IR photons 
or starlight. And the two-dimensional regression 
analysis further suggests L 11 energy density of the 
background optical/IR photon field and the stel- 
lar encounter rate/metallicity span a set of funda- 
mental planes (Hui et al. 2012). The unpulsed 
level inferred for M28 can be used to discrimi- 
nate which relation(s) can better predict the col- 
lective contribution. Using the best-fit parame- 
ters for these fundamental plane relations (Equa- 
tions 1-4 and Table 3 in Hui et al. 2011) and the 
updated GC parameters (Harris 1996; 2010 ver- 
sion), the relations involve metallicity and opti- 
cal/IR energy densities result in an estimate in 
a range of ~ (8 — 9) x 10 34 erg s _1 which is 
consistent with the observed off-pulse luminosity. 
On the other hand, the best-fit relations that in- 
volve the encounter rate result in an estimate of 
~ 2 x 10 35 erg s _1 which apparently overshoot the 
observed value. 

For further investigating this putative peri- 
odic signal, multi- wavelength observations are cer- 
tainly required. In particular, the phase-aligned 
X-ray/7— ray pulse profile will provide an impor- 
tant constraint for the high energy emission model. 
However, no existing X-ray timing data is avail- 
able for M28A in Fermi era. As the timing noise 
of M28A is quite strong in comparison with other 
MSPs and it possibly exhibited glitches, the phase- 
alignment of multi-wavelength light curves sub- 
jects to a lot of uncertainties. Therefore, follow- 
up timing observations in other wavelengths are 
encouraged for further investigations. 
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Tabic 1: Ephcmeris of PSR J1824-2452A adopted from Ray et al. (2008). 



Parameter 


Pulsar name 


J1824-2452A 


Right ascension, a 


18:24:32.00790550 


Declination, 5 


-24:52:10.8076448 


Pulse frequency, v (s _1 ) 


327.4056060517495439 


First derivative of pulse frequency, v (s~ 2 ) 


-1.735361869603xKT 13 


Epoch of frequency determination (MJD) 


53800 


Epoch of position determination (MJD) . . 


53800 


Solar system ephemeris model 


DE405 


Time system 


TDB 



5 



55800 




55600 



9, 55400- 
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Fig. 1. — Fermi LAT 7-ray weighted light curve (upper panel) and the phaseogram (lower panel) of M28A. A weight was 
assigned to each photon with the probability that it comes frogi M28 by using the task gtsrcprob in Fermi Science Tool. Two 
periods of rotation with a resolution of 40 phase bins per period is shown for clarity. The error bars of the light curve represent 
ltr Poisson uncertainties. The shaded regions define the on-pulse intervals for Peak 1 and Peak 2. 
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Table 2 

Phase-resolved spectroscopy of M28A 





Peak 1 


Peak 2 


off-pulse component 


PL fit 


r 


2.20 ±0.06 


2.18 ±0.09 


2.17 ±0.07 


/ph a 


(1.45 ±0.14) x 10~ 8 


(9.77 ± 1.20) x 10~ 9 


(1.02 ±0.13) x 10~ 8 


TS 


330 


303 


192 


PLE fit 


r 


0.68 ±0.32 


0.96 (fixed) 


1.26 ±0.27 


^cutoff 13 


1.16 ±0.27 


1.21 ±0.12 


2.05 ±0.71 


/ph a 


(1.01 ±0.14) x 10~ 8 


(8.29 ±0.79) x 10-° 


(8.42 ±1.24) x 10~ 9 


TS 


403 


368 


219 



a Photon flux in unit of photons cm 2 s 1 measured in the range of 0.2- 
300 GeV. 

a Cut-off energy in unit of GeV. 




Fig. 2. — Fermi LAT ph ase-resolved 7-ray count maps for events >0.2 GeV within 5° of the timing position of M28A 
(illustrated by the yellow cross). Top is north and left is east. The scale bar below shows the color scale of counts/pixel divided 
by the relevant phase interval. Left panel: Peak 1 region (i.e. 0.0< tj> <0.4). Middle panel: Peak 2 region (i.e. 0.55< <f> <0.75). 
Bight panel: Off-pulse region (i.e. 0.4< <j> <0.55 & 0.75< 4> <1). The point source locates in the northeast of these maps is 
2FGL Jf833. 6-2104 and the bright extended emission in the northwest is due to the diffuse 7— ray emission from the Galactic 
plane. 
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